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The surface dynamics and thermodynamics of metal nanowires are investigated in a continuum 
model. Competition between surface tension and electron-shell effects leads to a rich stability 
diagram, with fingers of stability extending to extremely high temperatures for certain magic con- 
ductance values. The linearized dynamics of the nanowire's surface are investigated, including both 
acoustic surface phonons and surface self-diffusion of atoms. On the stability boundary, the surface 
exhibits critical fluctuations, and the nanowire becomes inhomogeneous. Some stability fingers co- 
alesce at higher temperatures, or exhibit overhangs, leading to reentrant behavior. The nonlinear 
surface dynamics of unstable nanowires are also investigated in a single-mode approximation. We 
find evidence that some unstable nanowires do not break, but rather neck down to the next stable 
radius. 

PACS numbers: 68.35.Ja, 47.20.Dr, 61.46.+W, 68.65.La 



I. INTRODUCTION 

Metal wires play an essential role in all electrical cir- 
cuits, from power distribution between cities to inter- 
connects in integrated circuits. In today's technology, 
feature sizes down to approximately lOOnm are the state 
of the art, but current trends^ consistent with Moore's 
law^ extrapolate to lnm technology by 2020. A question 
of fundamental importance is whether metal will retain 
its role as the conductor of choice even at the ultimate 
limit of atomic-scale technology, or whether it must be 
replaced with more exotic conductors, such as carbon 
nanotubes. 3 

A macroscopic analysis of the mechanical properties 
of thin metal wires suggests that it might be difficult 
to fabricate wires thinner than a few thousand atoms 
in cross section: Consider a cylindrical wire of radius R 
and length L. The maximum tension that the wire can 
sustain before the onset of plastic flow is Fy — nR 2 cry , 
where cry is the yield strength. On the other hand, the 
force due to the surface tension cr s in a thin wire is 
F s = —nRo- s . If \F S \ > Fy, one would expect the wire 
to undergo plastic flow and, if L > 2irR, to break up 
under surface tension, as in the Rayleigh instability of a 
column of fluids This estimate gives a minimum radius 
for solidity, i? m i n = o~s/o~y. The parameters for several 
simple metals are given in Table [|] Plateau realized as 
early as 1873 that this surface-tension driven instability 
of a cylinder is unavoidable if cohesion is due solely to 
classical pairwise interactions between atoms£ 

A great deal of experimental evidence has accumu- 
lated over the past decade, however, indicating that 
metal nanowires considerably thinner than the above 
estimate can be fabricated by a number of different 
techniques ! 11 ' 12 : 13 ' 14 : 15 ' 16 : 17 ' 18 Even wires with lengths 
significantly exceeding their circumference were found to 
be remarkably stableji 3 ^ 1 ^ 7 " indicating that some new 



mechanism must intervene to prevent their breakup. 

An important technique which has been used to model 
the energetics of metal nanowires is classical molecular 
dynamics, 19 ' 20-21 ' 22-23,24 which utilizes short-ranged in- 
teratomic potentials optimized to fit the bulk properties 
of solids. This technique has had considerable success, 
including predicting the formation of metal nanocon- 
tacts in STM experiments 1 ^ and predicting novel, non- 
crystalline order in nanowires^ 2 " However, this approach, 
which neglects quantum-size effects, is unable to avoid 
the Rayleigh instability in long wiresi^2i2i24 

A clue to the resolution of this problem was provided 
by the observation of electron-shell structure in conduc- 
tance histograms of alkali metal nano contacts^ Like 
the surface tension, quantum-size effects arising from the 
confinement of the conduction electrons within the cross- 
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TABLE I: The yield strength cryj- surface energy a s £ and 
curvature energy 7.)— of various monovalent metals. The 
values2*"^ in the free-electron model, a s (fem) = ef^f/SOtt 
and 7 s (fem) = 4£f&f/457t 2 , are shown for comparison. For 
a wire of radius R < a s /ay, the stress due to surface tension 
exceeds ay, signalling a breakdown of macroscopic elasticity 
theory. The electrical conductance G m i n of a ballistic wire of 
radius R m i n = cr s /oy is shown in the rightmost column, in 
units of the conductance quantum Go = 2e 2 /h. Note that 
G/Gq is approximately equal to the number of atoms that fit 
within the cross section for monovalent metals. 
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FIG. 1: Stability of cylindrical metal nanowires as a function 
of radius and temperature. Shaded regions indicate stability 
with respect to small perturbations, A(Ro,T) > 0; unshaded 
regions denote unstable configurations, A(Ro,T) < 0. Here 
Tf is the Fermi temperature, &f the Fermi wavevector, and 
Ro the mean radius of the wire. The quantized conductance 
values of some of the stable wires are indicated. 



section of the wire become increasingly important as the 
wire is scaled down to atomic dimensions. In fact, a 
linear stability analysis 2 ^ of metal nanowires within the 
free-electron model found that the Rayleigh instability 
can be completely suppressed for certain favorable radii. 

In this article, we investigate the surface dynamics and 
thermodynamics of simple metal nanowires in a contin- 
uum approach, in order to shed further light on their 
unusual stability properties. The starting point for our 
analysis is the thermodynamic stability diagram shown 
in Fig. ^ Competition between surface tension and 
electron-shell effects leads to a complex landscape of sta- 
ble fingers and arches extending up to very high tem- 
peratures: wires whose electrical conductance is a magic 
number 1, 3, 6, 12, 17, 23,... times the conductance quan- 
tum G = 2e 2 /h are predicted to be stable with respect 
to small perturbations up to temperatures well above 
the bulk melting temperature Tm ~ .OlTp, where Tp 
is the Fermi temperature. This finding suggests that 
metal nanowires may be remarkably robust, which is 
cause for optimism about their potential for nanoelec- 
tronics applications. 26 Fig. ^ is akin to a phase diagram 
for metal nanowires; the nature of the different phases is 
revealed in this article through a study of the surface dy- 
namics for small perturbations about a cylinder. We find 
that the stable fingers correspond to homogeneous (i.e., 
translationally invariant) phases, while the intervening 
regions correspond to inhomogeneous phases. 



II. THE MODEL 

The continuum model we employ allows for an ana- 
lytical treatment of the long-wavelength surface modes 
used to characterize the different phases in Fig. ^ as 
well as a correct treatment of quantum-size effects, which 



are essential to stabilize long nanowires. The ionic de- 
grees of freedom of the wire are modeled as an incom- 
pressible, irrotational fluid, and the conduction electrons 
are treated as a Fermi gas confined within the wire by 
Dirichlet boundary conditions at the surface. Electron- 
electron interactions are included only at a macroscopic 
level (by requiring the wire to be electrically neutral), 
since it is well knowr&ifliSi that the leading mesoscopic 
shell-correction to the energy is independent of interac- 
tions. Calculations including interactions at the mean- 
field levelS&SS yield shell effects very similar to those in 
the free-electron modelASS 

Modelling the ionic degrees of freedom as a fluid is 
motivated by the argument presented in Table [fl which 
indicates that metal nanowires thinner than a few thou- 
sand atoms in cross section should be very plastic. The 
free-electron model for the conduction electrons is ap- 
propriate to describe electron-shell effects in monovalent 
metals, and as the values in Table ^ indicate, even de- 
scribes some macroscopic properties of alkali metals semi- 
quantitatively. Although the continuum approximation 
is not justified a priori in the limit of atomically-thin 
wires, this model is nonetheless justified a posteriori by 
its success in describing simple metal clusters^ of com- 
parable dimensions. Cohesion and quantum transport in 
gold nanocontacts were also successfully described with 
this modeli2fli22i2ii The directionality of bonding due to 
contributions from p-, d- and /-electrons is of course ab- 
sent from the free-electron model, as are element-specific 
effects, such as the tendency toward surface reconstruc- 
tion, which was argued to play an essential role in the 
formation of atomic chains^ Nonetheless, wires with 
G — Go are predicted to be very stable within the 
free-electron model, and the strength of a metallic bond 
in such a wire is significantly greater than that in the 

buik» 

An empirical justification for our continuum model 
comes from experimental results indicating that electron- 
shell effects dominate over ionic ordering in sufficiently 
thin alkali metali& and goldi& wires. Yanson et a/4 6 - 
found an interesting interplay between electron-shell ef- 
fects and atomic-shell effects in alkali metal nanocon- 
tacts. Electron-shell effects were found to be most im- 
portant in the lighter elements lithium and sodium, pre- 
sumably due to the larger Fermi energies of the conduc- 
tion electrons and the lighter, more mobile ions, while 
atomic-shell effects were most important in the heavier 
element potassium. A crossover from electron-shell struc- 
ture to atomic-shell structure in conductance histograms 
was found for conductance values G/Gq ~ 36 in potas- 
sium, while electron-shell effects were found to dominate 
even for G/Go > 100 in lithium. An intermediate be- 
havior was observed for sodium. Interestingly, the com- 
petition between the two effects was found to be history 
dependent. In a particular sequence of histograms ob- 
tained by cycling a potassium break junction, an evolu- 
tion from atomic-shell structure to electron-shell struc- 
ture was observed^ Most recently, a similar interplay 



between electron-shell structure and atomic-shell struc- 
ture was also observed in gold nanocontacts—Si 

These fascinating experimental results cry out for 
deeper theoretical investigations of the stability and 
structure of metal nanowires. While the geometry of 
the nanowires studied in Rcfs. was not directly 

determined, they may be rather short due to the fabri- 
cation method, so that the connection^ to the contacts 
may play an important role. In this article, we study the 
more theoretically tractable — and more technologically 
relevant — problem of the stability and surface dynamics 
of long metal nanowires. Our analysis sho uld be directly 
relevant for the nanowires studied in Refs. Il3ll7l 

III. LINEAR STABILITY ANALYSIS 



3 




(2,1) (3,1) (4,1) 




(4,2) (5,2) (6,2) 



FIG. 2: Classical periodic orbit*—*— of an electron in a plane 
perpendicular to the axis of the wire, labeled (v,w), where v 
is the number of vertices and w is the winding number. 



In our continuum model, the ionic degrees of free- 
dom are completely determined by the surface coordi- 
nates of the wire. Motivated by the fact that only modes 
which preserve axial symmetry participate in the surface- 
tension driven instability of a cylinder^ we restrict our 
consideration to axially-symmetric perturbations, 



R(z,t) = Ro 



(1) 



where R(z,t) is the radius of the wire at position z 
and time i, i?o is the unperturbed radius, and b(q,t) = 
b*(—q,t) are complex Fourier coefficients. Periodic 
boundary conditions are assumed for a wire of length 
L, so that q n = 2-kti/L, with n an integer bounded by 
|n| < N « fci? L I V (a lattice cutoff). Since the total num- 
ber of atoms comprising the nanowire is unchanged by 
the perturbation, 6(0, t) is related to the other b(q n , t) by 
volume conservation 



6(0, t) 



b 2 (0,t) 
2R Rq 



1 N 



(2) 



and may be eliminated. 

For small perturbations, the grand canonical potential 
of the electron gas is quadratic in the Fourier coefficients 
b(q, t), and determines the potential energy U of the ions 
in the Born-Oppcnhcimer approximation, 

TV 

U = Uo(R a ,T) + «(«»; Ro, T)\b(q n , t)\ 2 , (3) 



where Uq(Rq,T) is the potential energy of an unper- 
turbed cylinder, 



Uo(R,T) 
L 



itR 2 u + 2-kRus - tt 7s + V(R, T). (4) 



Here u is the macroscopic free energy density of the elec- 
tron gas, u s is the surface tension, 7 S is the surface curva- 
ture energy (c.f. Table|IJ), and V is a mesoscopic electron- 
shell correction. The mode stiffness a(q; Rq,T) has the 



following formS^ in the semiclassical approximation, valid 
for long-wavelength perturbations: 



a(q;R,T) 



-2ir<7 s /R + 2ir(a s R- ls )q 2 



\dli- RiUi 



8 2 1 B 

--5^)V(R,T), (5) 



where 



V(R,T) 



2s 



oo oo 



2 , > , — — cos(fc_FL l , to -3f7r/2). 



7T ' ' * ' V L v 

w—1 v—2w 



(6) 

The sum in Eq. JH| includes all classical periodic or- 
bits (v,w) in a disk billiard (see Fig. L vw = 
2vRsui(ttw/v) is the length of an orbit, the factor f vw 
I for v = 2w, 2 otherwise, accounts for the invari- 
ance under time-reversal symmetry of some orbits, and 
a vw {T) = T vw /siah.T vw (t vw = irk F L vw T/2T F ) is a 
temperature-dependent damping factor— 

Since d 2 a/dq 2 > for all physically meaningful radii, 
long-wavelength perturbations cost the least energy^ 
and the stability of the wire is determined by the sign 
of A(R ,T) = a(q = 0;R Ol T). For A(R ,T) > 0, a 
nanowire is stable with respect to all small perturba- 
tions, and is hence a metastable thermodynamic state. 
For A{Rq,T) < 0, the nanowire is unstable. The stabil- 
ity diagram so determined is shown in Fig. ^ In Fig. 
□ the values a s = a s (FEM) and 7 S = 7 S (FEM), ap- 
propriate for alkali metals, were used (c.f. Table 0. For 
larger values of a s (e.g. for noble metals), the maximum 
temperatures (in units of Tp) of the stable fingers are 
reduced somewhat, but the stability diagram is qualita- 
tively similar. 

Further insight into the stability criterion A > is 
provided by the identity 



A(R ,T) 



( d 2 



1 d \ U (Ro,T) 



\dR 2 



Ro dRr 



(7) 



The wire can lower its potential energy via a volume- 
conserving separation into thicker and thinner segments 



4 



if and only if A < 0. A < thus corresponds to an 
inhomogeneous phase, while A > corresponds to a ho- 
mogeneous phase. 

Our analysis of stability in terms of the convexity of 
the constrained energy functional is quite different from 
simply comparing the energy of cylinders of different 
radius ;SSiS2i2£ which does not address the fundamental 
question: whether any cylinder is stable. We also point 
out that for a sufficiently large system, the number of 
atoms is conserved — neglecting sublimation — and the de- 
pletion of atoms from a finite segment of wire^S can be 
described as a finite-wavelength perturbation of a larger 
system. 

Note that our stability analysis is carried out at fixed 
L. The tensile force necessary to fix the length of the wire 
is given by F — —dUo/dL [plus a small correction due 
to surface fluctuations, c.f. Eq. i|22|l ]. and was previously 
c alculat ed as a function of radius in this model in Refs. 
I^ lld30l Our stability analysis is thus appropriate to de- 
scribe nanowires under tensile stress, such as those stud- 
ied in the experiments of Refs. ITlTll2TITMT3IT5llT6lll7llT^ . 
The stability of a nanowire with free ends is an open 
question. 



IV. LINEARIZED SURFACE DYNAMICS 



where 1\ is the first-order modified Bessel function. The 
kinetic energy of the ionic medium is then given by 

K = Y J d 3 rV$*(r,t) -V$(r,t) 



JV 



= L^m(<7 n ,i? ) 



db(q n ,t) 



dt 



where pi is the ionic mass density, and 

2nRI {qR) 



m(q, R) = pi 



qh(qR) 



(12) 



(13) 



Details of the derivation of Eqs. (|12fl and (|13[) are given 
in Appendix [S] Combining Eqs. and 112f) yields a 
Hamiltonian for surface phonons, with frequencies 



j(q;R ,T) = 



a(g; Rp,T) 
m(q,R Q ) 



(14) 



Generically, uj(q) tx q as q — > due to the (/-dependence 
of m(q,Ro). Eq. (|14fl thus describes acoustic surface 
phonons. On the stability boundary ^4 = 0, one has 
uj(q) (x q 2 as q — > 0. For A < 0, lo is imaginary, and long- 
wavelength modes grow exponentially [see Fig.[5Ja)]. 



A. Surface phonons 



Surface diffusion 



We first consider inertial dynamics of the ionic back- 
ground. Assuming that the ionic medium is irrotational 
and incompressible^ its velocity distribution v(f, t) can 
be written in terms of a potential satisfying the Laplace 
equation 



V 2 $(r-,i) = 0, 



(8) 



where v(r,t) = — V$(r, t). The general solution to this 
equation with axial symmetry, which is regular at r = 0, 
can be written^ 

N 

*(f,t)=*(r,*,t) = d(q n ,t)I (q n r)e^ z , (9) 

n——N 

where Jo is the modified Bessel function of order zero and 
r is the distance of an ion from the z-axis. 

For small deformations, the relation between the co- 
efficients d(q n , t) in the expansion J^J and the Fourier 
coefficients b(q n ,t) of the surface perturbation J[J can be 
determined by the condition that the radial component 
of the velocity at the surface is 

Vr = —\r=R = (10) 

plus terms 0(b 3 ). Therefore, we have 

1 db(q n ,t) 



d(q n ,t) 



q n h{q n Ro) dt 



(11) 



The surface deformation also produces a gradient 
in chemical potential that drives the surface atoms to 
diffuse — a process likely to be important for large-scale 
deformations^ The surface current of atoms is given by 
Fick's law 



J 



PsDs 

k B T 



V/i, 



(15) 



where p s is the surface density of atoms and D s is the sur- 
face self-diffusion constant. Using the continuity equa- 
tion for the surface current, Eq. (|15|) can be converted 
into a (linearized) equation of motion for the profile 
R(z,t) 



dR(z,t) p s D s v a d 2 p 



dt 



k B T dz 2 ' 



(16) 



where v a = Sir 2 /hp is the volume of an atom. The chem- 
ical potential p of an atom is obtained by calculating the 
change in free energy with the addition of an atom at 
point Zo, 

p{z ,t) = U[R(z, t) + CS(z - z )] - U[R(z, t)], (17) 

where C = v a /2irR is chosen so that the volume of an 
atom is added. From Eq. J3J), one obtains 

p(z,t) = p + ^- y a(q n ;R ,T))b(q n ,t)e^ z , 



n— — N 



(18) 



5 



where /io(-Ro,T) is the chemical potential of the unper- 
turbed cylinder. Combining Eqs. (ftT))) and (fTHf> yields an 
equation of motion for the Fourier component b(q,t), 

^^l = -r(q;R Q ,T)b(q,t), (19) 

where the relaxation rate 

r(g; flo, T) = P ffi g 2 a(g; Ro ) T). (20) 

Thus, one finds that under surface diffusion alone, a per- 
turbed metastable wire relaxes exponentially toward a 
cylindrical shape. For a < 0, the mode grows exponen- 
tially. 
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FIG. 3: Root-mean-square fluctuations of the radius of 
metastable sodium nanowires of length L = 10J?o- 5R is 
undefined for unstable wires within the harmonic approxima- 
tion, and is not shown. 



C. Combined dynamics 

Combining inertial and diffusive processes, the lin- 
earized (classical) equation of motion for the surface 
modes is 



d 2 b(q,t) db(q,t) 



dt 2 



dt 



LO 2 (q)b(q,t)=0. (21) 



From the g-dependence of Eqs. I|14|) and l|20(l . one sees 
that T(q)/uj(q) — > as q — ► 0, indicating that diffusive 
processes can be neglected in this limit, at least for small 
deformations. In general, the relative time scales for in- 
ertial and diffusive dynamics depend on the value of D s . 
For this quasi-one-dimensional diffusion problem, one can 
estimate D s ~ {uj£,/p s ) exp(— E s /kgT), where loq is the 
Debye frequency and E s is the activation energy for sur- 
face diffusion, which is comparable to the energy of a 
single bond in the solid. With this estimate, one has 
u{q) 3> r(q) for all q, indicating that the surface phonons 
are underdamped. 

The picture of the surface dynamics of (meta)stable 
metal nanowires which emerges from this analysis is that 
there is a separation of timescales: on short timescales, 
the surface oscillates rapidly about the cylindrical equi- 
librium shape, while on much longer timescales, surface 
atoms diffuse irreversibly. 



V. CRITICAL SURFACE FLUCTUATIONS 

In the harmonic approximation [Eqs. JJJl and l|12l) ]. 

the total free energy £1(Rq,T) of the nanowire is given 
by the free energy of the unperturbed cylinder plus the 
Helmholtz free energy of the surface phonons, 

N 

n = U + J2 [ huJ n + 2k B T In (1 - e - f3huJ ")] , (22) 

n=l 

where uj n = Lo{q n \ R$,T) and (3 — l/fc^T. The equilib- 
rium tension in the wire is 



dL~ dL + pho 



(23) 



where the main c ontributi on —dUo/dL was previously 
calculated in Refs. I{"ll("l3(l and (JF'phonon is a small cor- 
rection that is singular at the stability boundaries, where 
the surface modes become soft. 

The softening of the surface modes on the stability 
boundaries leads to critical surface fluctuations. Given 
the stiffness (0 and frequency Ijl4|l of the surface modes, 
the mean-square thermal fluctuation SR. 2 of the radius of 
the nanowire can be calculated in the usual way^i 



1 N 

SR 2 ^((R-Ro) 2 ) = I Y / 



hw n [2f{Lj n ) + 1] 
a{q n ;R (h T) 



(24) 



where f(uj) = [exp(/3huj) — 1] _1 is the Planck distribu- 
tion. Note that, aside from a small quantum correc- 
tion, the magnitude of the surface fluctuations follows 
from the equipartition theorem applied to Eq. @, and 
is thus largely independent of the nature of the surface 
dynamics — whether inertial or diffusive. 

Fig. [3] shows SR for nanowires of finite length at room 
temperature as a function of their mean radius. The 
ionic mass and Fermi temperature were taken to be that 
of sodium. Within a metastable region, kpSR 1 
and is approximately independent of L, indicating that 
such wires are nearly atomically smooth at this temper- 
ature. The zero-point motion contributes roughly fifty 
percent of the surface fluctuation within a stable region 
for sodium nanowires at T = 300i^. With increasing 
temperature, the thermal contribution to 5R grows pro- 
portional to VT, according to the equipartition theorem. 
Note that the harmonic approximation is expected to 
break down when kpSR ~ 1. 

The surface fluctuation SR exhibits sharp peaks at the 
stability boundaries, reaching the value 



SR 2 



k B T 



A(R ,T)=o 47T 3 (ai? - 7 ) 



L. 



(25) 



plus a small quantum correction. SR thus scales as L 1 / 2 
on the stability boundary, like the finite-size scaling at 
the roughening transition of a planar interfaced At T = 
0, SR remains small and approximately independent of L 



6 



on the stability boundary. The absence of critical surface 
fluctuations at T = is also consistent with the behavior 
of planar interfaces^ 

The stability boundary A(Rq,T) — defines a 
(multiple- valued) critical temperature T c = T c (Rq) as 
a function of the mean radius, or alternatively a criti- 
cal mean radius R c — R C (T) as a function of tempera- 
ture (see Fig. Within the harmonic approximation, 
5R grows with an exponent v = -1/4 as i?o -> i? c or 
T — > T c , as expected from the Ornstein-Zernicke fluctua- 
tion theory. This critical behavior, which is cut off when 
5R approaches the value given in Eq. I)25[l. is illustrated 
in Fig. El 

One can also study the time-dependence 4 ^ of 6R for an 
initially cylindrical wire undergoing thermal fluctuations. 
Using the classical equipartition theorem, for each mode 
we have 



2ksT sin 2 (uj n t) 
L a(q n ;R ,T) 



(26) 



where sm 2 (co n t) describes the standing capillary waves. 
The surface fluctuations then grow as a function of time 
according to 



SR 2 (t) 



Ak B T 



N 

E 



sin 2 (oj n t) 
J a(q n ;R ,T)' 



(27) 



On the stability boundaries A(Rq, T) = 0, the dispersion 
relation to(q) ~ q 2 , and one finds asymptotically SR(t) ~ 
t 1 / 4 for times ui^ 1 t uj^ 1 . The dynamic exponent 
2 = 1/4 differs from that of a planar interface 4 ^ due to 
the different dispersion relation for the surface modes. 

However, all these scaling relations hold only in a lim- 
ited range, since the asymptotic limit SR — > oo charac- 
terizing the roughening transitions 41 *^ is unphysical in 
nanowires, due to their finite radius. On the stability 
boundary, the surface does not roughen in a thermody- 
namic sense, but the nanowire does become inhomoge- 
neous. 



VI. REENTRANT BEHAVIOR 

Perhaps most interesting is the reentrant behavior oc- 
curing on the arches and overhangs in the stability di- 
agram, Fig. ^ For instance, a wire with kpR^ — 19 is 
metastable and homogeneous in the temperature interval 
T cl <T < T c2 , with T cl w .0072T F and T c2 w .046T F . 
The surface exhibits critical fluctuations as T — > or 
T — > T C 2, at which points the wire makes a transition to 
an inhomogeneous phase. The transition at T c2 is con- 
ventional, in the sense that the inhomogeneous phase is 
the high-entropy phase. However, the inhomogeneous 
phase below T c \ has lower entropy than the homoge- 
neous phase above T c \. Fig.^Ja) shows the total entropy 
S = —dfl/dT of the nanowire as a function of tempera- 
ture, including both electron and phonon contributions, 
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FIG. 4: (a) The total entropy per unit length of a nanowire 
with kpRt, — 19 versus temperature. The ionic mass was 
taken to be that of sodium, (b) The electron-shell potential 
V(Ro,T), from Eq. ©. 



where fi is given by Eq. I|22|) . The electronic entropy is 
regular at the critical points, but the phonon entropy is 
singular in the harmonic approximation, due to the emer- 
gence of soft surface modes. The singular contribution to 
the phonon entropy is 



■£/(<■ 

n=l 



duj{q n ; Rp,T) 
dT 



(28) 



indicating that the softening of a phonon mode with de- 
creasing temperature indeed leads to a decrease in en- 
tropy. 

To understand the counterintuitive behavior at T c i, it 
is useful to consider the electron-shell correction V to 
the energy of the wire, shown in Fig. E[b). Above T c \, 
V has a single broad minimum near kpRo = 19, but as 
the temperature is lowered, and the fine structure in the 
shell potential emerges, this single minimum splits into 
two minima at kpRo — 18.75 and kpRo = 19.2. To lower 
its free energy, the system would like to fall into one of 
these two minima, but due to volume conservation, such 
a global change is not possible. The wire thus undergoes 
phase separation into thick and thin segments^ 



VII. UNSTABLE WIRES 

Finally, let us discuss the dynamics of unstable wires. 
Fig. a) shows the real and imaginary parts of the sur- 
face phonon frequency versus wavevector for a typical 
unstable wire. One mode b(q m ,t) grows exponentially 
faster than all others in the harmonic approximation, and 
thus may be expected to dominate. For a single Fourier 
component b(q m ), the potential energy U[b(q m )] of the 
nanowire can be evaluated for arbitrarily large b using 
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semiclassical perturbation theory. U may be expanded 
semiclassically a.s££l 

U = uV + a s S - 7s C + SU. (29) 

The volume V, surface area 5, and integrated mean cur- 
vature C of the nanowire can be calculated for arbitrary 



5U[b(q m )] _ 2ef_ a vw (T)f vw 

L ~ it ^ ^ v 2 L vw 

w—l v—Zw 



where 8 VW = k F L vw (l + b(0)/Ro) ~ 3vir/2, <p vw = 
2kFL vw b(q m )/ Ro, Ji is the i-th order Bessel function, 
and b(0) is related to b{q m ) by Eq. @. The result is 
shown in Fig. GJb). Although the wire is unstable to 
breakup under a hypothetical long-wavelength perturba- 
tion, the energy of the fastest growing mode reaches a 
minimum at a finite amplitude, suggesting that the sur- 
face deformation saturates, and that the wire does not 
break up, but rather necks down to the next stable ra- 
dius. A similar scenario is predicted under the diffusive 
dynamics of Eq. 1)190. An explicit nonlinear dynamical 
simulation 35 confirms these predictions. 
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FIG. 5: (a) The dispersion relation for the surface modes of 
an unstable nanowire with fcp-Ro = 8.95 at T = 0. The ionic 
mass was taken to be that of sodium, (b) The potential en- 
ergy U[b(q)] of the same nanowire for sinusoidal deformations 
with q = and q = q m — 0.208/cf, respectively. Note that the 
energy Uo of a straight wire has been subtracted. The differ- 
ent energies of these two modes is mainly due to the increased 
surface energy at finite q. 



deformations by simple geometric considerations. Us- 
ing semiclassical perturbation theory^Mil 5 . the electron- 
shell correction SU can again be expressed in terms of the 
classical periodic orbits of a disk billiard, leading to an 
expression similar to Eq. JfjJ): 



(30) 

I 

VIII. CONCLUSIONS 

The stability and surface dynamics of metal nanowires 
were investigated in a continuum approach, including 
electron-shell effects. A thermodynamic phase diagram 
for jellium nanowires was derived, which predicts that 
cylindrical wires with certain "magic" conductance val- 
ues are stable with respect to small perturbations up to 
remarkably high temperatures. On the stability bound- 
ary, the surface exhibits critical fluctuations, and the 
nanowire becomes inhomogeneous. Both surface phonons 
and surface self-diffusion of atoms were included in the 
linearized surface dynamics. It was found that in- 
ertial dynamics (phonons) always dominate the long- 
wavelength behavior, including the critical points. (It 
must be emphasized, however, that this conclusion holds 
only for small perturbations of the surface. The irre- 
versible diffusion of surface atoms is undoubtedly crucial 
for large-scale surface deformations £ty A novel reentrant 
behavior was predicted, in which a straight wire is stabi- 
lized at intermediate temperatures, but undergoes phase 
separation into thick and thin segments as the temper- 
ature is lowered. Finally, for unstable wires, the surface 
deformation was found to grow exponentially, dominated 
by a single Fourier component, and to saturate at a finite 
amplitude, indicating that unstable wires may not break, 
but rather neck down to the next stable radius. 

The results presented in this article should be di- 
rectly relevant for nanowires made of monovalent met- 
als, especially the alkali metals and gold, for which 
there is clear experimental evidence of electron-shell 
cffectsii£ii§ii& Moreover, this simple model may provide 
qualitative insight into the generic surface properties 
of metal nanowires, which could guide investigations of 
more realistic, material-specific models. 
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APPENDIX A: THE IONIC KINETIC ENERGY 

Here we present some details of the derivation of the 
expression (I12fl for the kinetic energy of the ions. Insert- 
ing Eq. into the first line of Eq. (|12(l and performing 
the z-integral, one obtains 



i-Ro 

K = tt Pi L V / rdr [qllg(q n r) + qlrf(q n r)] 
„ Jo 



\d(q n ,t)\ 2 . 
(Al) 



dx 



[xl m (x)l' m {x)\ = X 



(A2) 



the radial integral in Eq. (|A1|) may be performed, leading 
to the result 



K = v Pl L qnRoIo(qnRo)h{q n Ro)\d(q n )\ 2 ■ (A3) 



n=-N 



Finally, eliminating d(q n ,t) from Eq. (|A3|) using the re- 
lation , one obtains the second line of Eq. (11211 . 
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